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Abstract. The blazar lES 1959-1-650 was observed twice by BeppoSAX in September 2001 simultaneously with 
optical observations. We report here the X-ray data together with the optical, Rc magnitude, light curve since 
August 1995. The BeppoSAX observations were triggered by an active X-ray status of the source. The X-ray 
spectra are brighter than the previously published X-ray observations, although the source was in an even higher 
state a few months later, as monitored by the ASM onboard RossiXTE, when it was also detected to flare in the 
TeV band. Our X-ray spectra are well represented by a continuosly curved model up to 45 keV and are interpreted 
as synchrotron emission, with the peak moving to higher energies. This is also confirmed by the slope of the X-ray 
spectrum which is harder than in previous observations. Based on our optical and X-ray data, the synchrotron 
peak turns out to be in the range 0.1-0.7 keV. We compare our data with non simultaneous radio to TeV data and 
model the spectral energy distribution with a homogeneous, one-zone synchrotron inverse Compton model. We 
derive physical parameters that are typical of low power High Energy peaked Blazar, characterised by a relatively 
large beaming factor, low luminosity and absence of external seed photons. 



Key words, galaxies: BL Lacertae objects: general - galax- 
ies: BL Lacertae objects: individual: lES 1959-1-650 - X- 
Rays: galaxies 

1. Introduction 

Blazars constitute the most extreme class of active galactic 
nuclei: they are radio loud AGN characterized by high lu- 
minosity, with rapid and high amplitude variability. They 
are usually divided into Flat Spectrum Radio Quasars and 
BL Lac objects, according to the presence or not of broad 
spectral features. Another blazar distinguishing property 
is the spectral energy distribution which is marked by two 
broad peaks (von Montigny et al. 1995): the first, extend- 
ing from the radio to the UV/X-ray band is usually inter- 
preted as synchrotron emission, while the second, which 
reaches the 7-ray band (sometimes TeV frequencies) is ex- 
plained as an inverse Compton component. Based on this, 
BL Lacs are further split into two subclasses distinguished 
by the synchrotron peak frequency: LBL and HBL, Low 
and High Energy peaked BL Lacs, respectively (Padovani 
& Giommi 1995). Various models have been proposed to 
explain the, still debated, origin of the comptonized pho- 
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tons: they could be the synchrotron photons themselves 
(Maraschi, Ghisellini & Celotti 1992), maybe partly repro- 
cessed by the broad line region (Ghisellini & Madau 1996) 
or external ambient photons emitted by the accretion disk 
(Dermer & Shhckeiser 1993) possibly reprocessed to some 
extent by the broad line region (Sikora, Begelmann & Rees 
1994) or by a dusty torus (Blazejowski et al. 2000). 



A wide spectral coverage is therefore necessary to un- 
derstand the physics of these objects: accurate simulta- 
neous multiwavelength spectra and light curves allow us 
to constrain the mechanism at work and the geometry 
of the emitting region. HBL in particular have gained 
much attention in the last decade since many of them 
have proved to be TeV emitters: the strongest ones are 
the thoroughly investigated Mkn421 (Punch et al. 1992) 
and MknSOl (Quinn et al. 1996). Another three sources 
have then been detected by ground-based Cherenkov tele- 
scopes: lES 2344-^514 (Catanese et al. 1998), PKS2155- 
304 (Chadwick et al. 1999) and H 1426-^428 (Horan et al. 
2002). On the basis of X-ray and radio data, Costamante 
& GhiseUini (2002) selected a catalogue of BL Lac objects 
which are most likely to be TeV emitters. One of the most 
promising objects in this catalogue is lES 1959-1-650, as 
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suggested also by Stecker et al. (1996). Indeed, this source 
has now been detected in the TeV band (see below). 

lES 1959+650 was discovered in the radio band as 
part of a 4.85 GHz survey performed with the 91 m 
NRAO Green Bank telescope (Gregory & Condon 1991; 
Becker, White & Edwards 1991). Subsequently the source 
was observed also in the optical band where it displayed 
large and fast flux variations (e.g. Villata et al. 2000). 
It shows a complex structure composed by an elliptical 
galaxy (M^ = —23) plus a disc and an absorption dust 
lane (Heidt et al. 1999). Its redshift (z=0.048) was derived 
by Perlman et al. (1996) from a spectrum obtained at the 
2.1m telescope at Kitt Peak. 

The first X-ray measurement of 1959+650 was per- 
formed by Einstein-IPC during the Slew Survey (Elvis 
et al. 1992). Subsequently, the source was observed by 
ROSAT in 1996, by BeppoSAX in 1997 (Beckmann 2000; 
Beckmann et al. 2002) and finally by USA and RXTE 
during 2000 (Giebels et al. 2002). On the basis of a X- 
ray/radio versus X-ray/optical color-color diagram the 
source was classified as a BL Lac object by Schachter et al. 
(1993). According to the new schemes, the source belongs 
to the HBL class. 

In the 7-ray band, it was observed in 1995 by EGRET 
on board the CGRO and by the Whipple observatory, but 
the source was not detected and only upper limits could 
be set (Weekes et al. 1996); subsequent EGRET obser- 
vations allowed Hartmann et al. (1999) to put it in the 
third EGRET catalog, with an average flux of 1.8 x lO^'' 
photons cm^'^ s^"'^ at energies above 100 MeV. 

During the spring and the summer of 1998, it was ob- 
served intensively with the Utah Seven Telescope Array 
detector which recorded a 3.9cr total significance above 
600 GeV. However, during two epochs, the source was seen 
above 5a significance (Nishiyama et al. 1999). These data 
have been confirmed by the HEGRA and the Whipple 
teams. After 94 hours on-source collected in 2000 and 
2001, the HEGRA team detected the source at 8% of the 
Crab fiux with a 5.4 a significance, while from May to July 
2002 the source was detected in a much brighter state, up 
to and above the Crab flux (Horns et al. 2002, Aharonian 
et al. 2003). Also the Whipple team observed the source 
for 39 hours between May and July 2002, reporting 
a mean flux of 0.64 ± 0.03 Crab above 600 GeV, with a 
significance of 20 a (Holder et al. 2003a, b). Thus, during 
the May- July 2002 HEGRA and Whipple observations, 
lES 1959+650 was one of the strongest TeV sources in 
the sky. 

lES 1959+650 is therefore one of the most interest- 
ing and frequently observed sources of recent years. In 
this paper we present two BeppoSAK ToO observations 
as part of a project for observing active state blazars as 
detected in optical. X-ray or TeV bands. In particular, 
these observations were triggered by the active X-ray state 
of lES 1959+650 as observed by the RossiXTE All Sky 
Monitor (ASM) and are the final observations of this suc- 
cessful BeppoSAX project. 



2. Observations and data reduction 

The Italian-Dutch satellite BeppoSAX has proved par- 
ticularly useful for studying blazars because of its ex- 
tremely large spectral range (0.1-200 keV) which allowed 
the detection of the transition between synchrotron and 
Compton emission (Tagliaferri et al. 2000; Giommi et al. 
2000; Ravasio et al. 2002) and the interesting comparison 
of simultaneous soft-X versus hard-X behaviour (Fossati 
et al. 2000a; Zhang et al. 2002). 

For a detailed description of the mission we refer to 
Boella et al. (1997 and references therein). In our work 
we have analyzed the data of different co-aligned Narrow 
Field Instruments aboard the satellite: the Low Energy 
Concentrator Spectrometer (LEGS; 0.1-10 keV), the two 
(originally three) identical Medium Energy Concentrator 
Spectrometers (MFCS; 1.5-10 keV) still operating and the 
passively collimated Phoswich Detector System (PDS; 13- 
300 keV). 

Since lES 1959+650 was observed to be in a high X- 
ray state by the iZossiXTE-ASM during September 2001, 
we triggered a BeppoSAX ToO observation, scheduled to 
last 5 X lO'' sec. However, because of technical problems, 
it was stopped after ~ 5 x 10'^ s. Thus, we were given a 
second opportunity and some days later a new pointing 
was performed. In Tabled we show the log of BeppoSAX 
observations, together with exposures and count rates for 
each instrument. 

We based our analysis on linearized and cleaned event 
files available from the online archive of the BeppoSAX 
Science Data Center (Giommi & Fiore 1998): the events 
from the two operating MECS are merged together to im- 
prove the photon statistic. Using XSELECT V2.0, we ex- 
tracted light curves and spectra from the event files se- 
lecting circular regions centered on the source of 8 and 
4 arcmin radii for LEGS and MECS respectively, as sug- 
gested in Fiore et al. (1999). We also extracted events from 
off-source regions to test the background stability. For 
the spectral analysis we preferred to use background files 
accumulated from long blank field exposures and avail- 
able from the SDC public ftp site: this choice is motivated 
by the non uniformity of LEGS and MECS backgrounds 
across the detectors (Fiore et al. 1999; Parmar et al. 1999). 
In this way the source and background spectra are ex- 
tracted from the same detector area. For the spectral and 
temporal analysis we used XSPEC 11.0.1 and XRONOS 
5.18 packages, respectively. 



3. Spectral analysis 

In order to reproduce the whole X-ray spectrum we fit- 
ted the data of the LEGS, MECS and PDS together: to 
account for the uncertainties in the intercalibration of the 
instruments it is necessary to use constant rescaling fac- 
tors. The LECS/MECS normalization acceptable range is 
[0.67-1], while the PDS/MECS one is [0.77-0.93] (Fiore et 
al. 1999). In order to reduce the number of fitting param- 
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LEGS 


MEGS 


PDS 


Date 


exposure count rate^ 
(s) 


exposure count rate'' 
(s) 


exposure count rate"^ 
(s) 


25 September 2001 


5.8x10'' 1.12 ±0.015 


7.2x10^ 1.43 ±0.016 


2.97x10^ 0.23 ±0.099 


28-29 September 2001 


2.6x10* 1.23 ±0.008 


4.8x10** 1.75 ±0.007 


2.3 x 10"* 0.35 ±0.036 



Table 1. Log of BeppoSAX observations. ^ 0.1-10 keV; ^ 1.5-10 keV; " 15-100 keV. 



eters and to be consistent with previous work, we fixed 
them at 0.72 and 0.9, respectively. 

3.1. 25 September 2001 

Although the first observation was very short the source 
was detected even by the PDS up to 35 keV: the good 
channels were then grouped in 3 bins. We first analyzed 
LEGS + MEGS spectra and then we added the PDS 
data: since for this observation and for that of the 28- 
29 September, the addition of PDS data does not change 
significantly the best fit parameters, in our discussion we 
will always refer to the analysis performed on the full set 
of LEGS-MECS-PDS spectra. 

We fitted the data with a power-law or a broken 
power-law model, first fixing the absorption parameter 
Nh to the Galactic value (Nh = 1.0 x 10^^ cm~^ as deter- 
mined from 21 cm measurements by Dickey & Lockman 
(1990)) and then leaving it free: the best fit parameters 
of each model are given in Table [21 together with the flux 
at 1 keV and the integrated [2-10] keV flux. Using the 
Galactic absorption value we always obtained bad fits: 
large negative residuals toward low frequencies are ev- 
ident. With the Nh free to vary, instead, we obtained 
good fits to the spectrum both with a simple power-law 
( energy index a = 1.45 ± 0.05, Nr = 2.1 x 10^^ cm-^; 
Xr/d-o.f. — 0.7/32) and a steepening broken power-law 
model (Nh = 1.8 x lO^i cm-^; xlld.o.f. = 0.6/30). The 
best-fit parameters of the broken power-law model are 
ai = 1.3 ± 0.1 and a2 = 1.5 ± 0.1 with the break located 
at 2.7 keV. The values do not allow us to distinguish 
between these spectral models which are both acceptable. 
Note, however, that in both cases the resulting Nh values 
are significantly higher than the Galactic value. 

These results are consistent with previous X-ray ob- 
servations of the source. Beckmann et al. (2002) an- 
alyzed the 1997 BeppoSAX and the ROSAT AU Sky 
Survey data finding a Nh value in excess of the Galactic 
value: more precisely, using a power-law model they found 
Nh = 2.55 X 10^^ cm~^ for the BeppoSAX spectrum and 
Nh = 1.6 X 10^1 cm~2 for the i?05^T spectrum. For the 
BeppoSAX observation they also fitted a broken power- 
law model finding an intermediate value. This extra ab- 
sorption required to fit the X-ray data could be due ei- 
ther to a real intrinsic absorption at the source or to an 
artifact due to a progressive steepening of the spectrum 
that can not be reproduced by a simple power-law or bro- 
ken power-law model. The assumption of an absorption 



higher than the Galactic value would not be necessary if 
a model which can intrinsically account for a progressive 
steepening of the spectrum is used. Therefore, we tried to 
reproduce the BeppoSAX spectrum with the continuously 
curved model described in Fossati et al. (2000b): 

FiE) = KE-"— [1 + (— (1) 
Eb 

where a_oo and a+oo are the asymptotic values of the en- 
ergy indices, while Eb and / are the parameters that de- 
termine the spectral bending. We applied this model to the 
whole LEGS-MECS-PDS spectrum fixing the absorption 
parameter to the Galactic value. First we found the value 
of the parameter / that minimizes the when all the 
other parameters are left free to vary obtainiong / = 1.1. 
Then, we repeated the fitting a few times to derive the 
spectral indices at several reference energies (0.5, 1, 10, 35 
keV) , checking that each time the same values for the vary- 
ing parameters were obtained (see Fossati et al. 2000b for 
more details of this procedure). The data are well fitted, 
with only the Galactic absorption, {xl/d.o.f. = 0.6/31) 
by a curve which is flat below 1 keV and steepens to- 
ward higher energies; the model reaches a slope a = 1 at 
1.2 ± 0.1 keV. In TableOwe give the best-fit parameters 
of the model. 

A possible inconvenience of this model is that it is so 
fiexible to adjust the parameters for a wide rage of spec- 
tral curvatures and then it results poorly sensitive to the 
actual Nh value. Therefore, we also considered another 
curved model using a logarithmic parabola, which pro- 
vides a reasonable representation of the wide band spectral 
distribution for the synchrotron component of blazars: 

F{E) = K {E/Eiy^''+^ Log(E/E,)) ^2) 

This model has recently been used to fit the optical- X-ray 
spectra of various objects (Giommi et al. 2002; Massaro et 
al. 2003a, b) . In our computations we fixed the Ex value at 
1 keV. With the Nh fixed to the Galactic value, this model 
provides a good fit to the data although not as good as 
the previous curved model. A lower x^ is obtained if we 
leave the Nh free to vary, the resulting best-fit value is 
very similar to the one obtained with the broken power- 
law model. A summary of the best-fit parameters is given 
in Table 13. 

3.2. 28-29 September 2001 

The second pointing had an exposure about a factor of 6 
longer than the previous one. Furthermore lES 1959+650 
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power law broken power law 




1 10 1 10 

channel energy (keV) channel energy (keV) 



Fig. 1. lES 1959+650 BeppoSAX spectrum of 28-29 September 2001. Left panel: LEGS + MEGS + PDS spectrum 
fitted with a power-law plus intrinsic absorption model: positive residuals are evident below 0.8 keV, while there are 
negative residuals in the PDS range. Right panel: LEGS + MEGS + PDS spectrum fitted with a broken power-law 
plus intrinsic absorption model provides a better fit, although the PDS data still lie below the model. 



was ^^20% brighter than in the first observation (see 
Table 1). We therefore had data with much higher counts, 
useful for a more accurate investigation of the spectral 
shape. The PDS was able to detect the source up to 
slightly higher frequencies (~ 45 keV) than the previ- 
ous observation. We again grouped the PDS data in 3 
bins. Also in this case we will refer to the analysis per- 
formed on LEGS+MEGS+PDS data as a whole, repeat- 
ing the fitting process with the absorption parameter 
fixed to the Galactic value or free to vary. Spectral fits 
using power-law or broken power-law models with the 
Galactic Nh are not acceptable because of the rather 
high values (see Table . With a free Nh we again 
find values that exceed the galactic value and are con- 
sistent with those of the September 25 observation and 
of Beckmann et al. (2002). However, in this case the 
BeppoSAX spectrum is also badly fitted by a single power- 
law model (Xr/d.o.f.=2.3/32) with the Nh free to vary. 
As can be clearly seen in Fig. ^ this model leaves pos- 
itive residuals toward low energies and negative ones in 
the PDS range: the spectrum seems to be continuously 
steepening. Leaving the absorption parameter free, a bro- 
ken power-law model with a break at E;, ~ 2.6 keV 
(Aa ^ 0.26) provides a better representation to the 
data (xr/d.o.f.=1.0/30), although the negative residuals 
at higher energies are still not completely removed (Fig. 

Similarly to what was obtained in the analysis of the 
previous observation, to fit BeppoSAX data with a bro- 
ken power-law model we need a column density higher 
than the Galactic value. We then fitted our data using 
the curved model of Fossati et al. (2000b). The data are 
well fitted with a curvature parameter / — 1.25 and a 
fixed Galactic absorption value {x^/d.o.f. ~ 1.1/31). In 
this case the negative residuals left in the PDS range by 
the broken power-law model disappear (see Fig. [21). Best- 
fit parameters are given in Table |31 Similarly to what was 



observed some days before, the model is hard below 1 keV 
and steepens toward higher energies up to a^^keV = 1-48; 
the turning point where the model reaches a = 1 is shifted 
to a slightly larger frequency (1.55 ± 0.05 keV). Finally, 
we applied the parabolic model, but this time a good fit 
could be obtained only with the Nh free to vary, finding an 
absorption value similar to that obtained with the broken 
power-law model (see Table . 

To summarize our results, on September 2001, 
lES 1959+650 X-ray spectra were convex and could be 
fitted either with a steepening broken power-law or a 
parabolic model with an interstellar absorption exceed- 
ing the Galactic value. Taking into account that optical 
images of the host galaxy (see Sect. 5) show that it is 
crossed by a strong dust lane, it is well possible that this 
additional Nh is due to some contribution within the host 
galaxy. The higher value of Nh derived from the X-ray 
data is in line with the E(B-V) values derived from optical 
observations (see section 5). The only model that fits the 
data with the Nh value fixed at the Galactic value is the 
curved model proposed by Fossati et al. (2000b). However, 
the X-ray spectral shape that one derives in this case is 
not very plausible in the overall SED of the source, with 
a synchrotron peak above 1 keV. Thus, we think that this 
representation, with the Nh fixed at the Galactic value, 
is not realistic. To conclude, our X-ray data imply an ab- 
sorption higher than the Galactic value and a peak of the 
X-ray emission in the 0.1-0.7 keV range. 

3.3. Archival X-ray data 

It is interesting to compare the 2001 BeppoSAX spec- 
tra to those previously measured by other X-ray mis- 
sions: Einstein, ROSAT (1996), BeppoSAX itself (1997), 
RXTE and USA (2000). The results of these observa- 
tions are summarised in Tabled During Einstein, ROSAT 
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model 



Nh 

l22 



(xlO^^ cm-^) 



keV 



a2 



FlkeV 



25 September 2001 



LEGS + MEGS + PDS 



X2r-/d.0.f. 



power-law 

power-law 
broken power law 
broken power law 



0.1 

0.21 ±0.03 
0.1 

0.18 ±0.0;.! 



1.25 

1.45 ± 0.05 

± 0.1 



1.9 



+ 0.5 

-().f< 



1.45 ±0.1 
1..") ± 0.1 



31.3 

40.3 
29.5 
37.2 



28-29 September 2001 



LEGS + MEGS + PDS 



8.6 

8.3 

8.3 
8.3 



4.95/33 

0.70/32 
1.20/31 
0.55/30 



power-law 
power law 
broken power-law 
broken power-law 


0.1 

0.22 

0.1 

n 1C-003 
0-lD_,.o.o2 


1.15 
1.3 

0.85 2.1 
1.1 ±0.1 2.6t°l 


1.3 

1 oc+0.03 
J--oD_0.04 


33.5 
42.1 
30.5 
37.0 


10.6 
10.4 
10.4 
10.4 


24/33 
2.3/32 
2.8/31 
1.0/30 


Table 2. Best fit model parameters. * F2- 


■lOkeV flux in unit of 10~ 


"^^ erg cm" 


-2 s-i. 







LEGS + MEGS + PDS: curved model 


25 September 2001 










^ 2-U)keV 


Xr/(1.0.f 




0.1) ± 0.1 


'-■''''-il.Dr, 


^■('-(),ur, 


8.:i 


0.()/:il 


28-29 September 2001 


0.2 ±0.2 


0.74 ± 0.03 


1.43 ± 0.03 


i-48+;;:;;:i 7.8 


10.4 


1.1/31 




LEGS 


± MEGS ± 


I'DS: parabolic iiiiulol 










25 September 2001 






Nh 


a 


b 


FlfceV F'o5_2keV 








'^2-lOkeV 


Xr/d.O.f. 




















0.17 ±0.06 


1.24 ± 0.27 


0.18 ±0.25 


36.9 7.49 


8.26 


0.63/31 


FIXED 


1.925 ± 0.075 


0.465 ±0.1 


30.5 7.43 


8.22 


0.96/32 


28-29 September 2001 


O.K) ^ ().():! 


1.0() ± 0.10 


0.21 = 0.09 


37.0 7.82 


10.1 


1.18/31 


FIXED 


1.83 


0.4 


31.7 7.82 


10.4 


2.4/32 



Table 3. Best fit parameters for the curved and parabolic models of LECS-MECS-PDS spectrum. * fluxes are 

in units of 10^^^ erg cm^^ s^"'^. For the curved model, the absorption parameter is fixed to the Galactic value of 
Nh = 1-0 X 10^^ cm~2; while the best-fit curvature parameters are / = 1.1 and / = 1.25 for the first and the second 
observation, respectively, reaching a slope a = 1 at 1.2 ± 0.1 and 1.55 ± 0.05 keV for the two spectra. 



and the 1997 BeppoSAX observations, the source had 
a flux of F2-iQkeV ~ 10~^^ erg cm~^ s~^. In particu- 
lar, Einstein measured a 2 keV flux which was ^ 40% 
and ~ 60% higher than the 1997 BeppoSAX and 1996 
ROSAT measurements, respectively. lES 1959+650 was 
monitored by RXTE and by USA through the summer 
and the autumn of 2000. The source was in a higher state 
(F2_iofcey ~ lO"-^*^ erg cm"^ s~^), with a flare reaching a 
2-10 keV flux of 2.3 x 10"^° erg cm^^ g-i (^4 November 
2000, USA; Giebels et al. 2002). The source was in a 
high state also during our observations, similar to those 
measured by RXTE and a little fainter than the flare of 
November 2000, monitored by USA. 

These X-ray spectra are always soft (a > 1) and can be 
explained as synchrotron emission. The BeppoSAX spec- 
tra of 2001 arc the hardest measured so far (below 3 kcV), 
as if the synchrotron component were peaking at higher 



energies with respect to the other historical observations 
(see Discussion). 



4. Temporal analysis 

We examined the observation of September 25 for short 
time variability, rebinning the light curves with a resolu- 
tion of 600 seconds and selecting only bin intervals with an 
effective exposure greater than 20% (i.e. all bins that have 
data for less than 20bin are rejected). Because the spec- 
trum seems to steepen at ~ 3 keV, we analysed the LEGS 
and MEGS light curves in two different energy ranges: 
[0.3-3] keV LEGS and [3-10] keV MEGS bands. During 
this short period, the observation was interrupted after 
^ 7000 s, the source seemed to be quiet: both LEGS [0.3- 
3] keV and MEGS [3-10] keV light curves are well fitted 
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Instrument 


Date 


a 


F2-10fceV 

(10"" erg cm"^ s"^) 


ROSAT 


31 March 1996 


J^- '"-0.21 


1.04 


BeppoSAX 


4 May 1997 


1.64 ±0.08 


1.29 


RXTE 


July-September 2000 


1.38-1.68 


8.4-14 


USA 


October-November 2000 




4.87-22.6 


BeppoSAX 


25 September 2001 


1.31" 


8.29 


BeppoSAX 


28-29 September 2001 


1.12" 


10.42 



Table 4. Historical X-ray observations of lES 1959-1-650. " soft X-ray spectral index. 



curved model 



36-39/9/2001 



28-29 September 2001 



channel energy (keV) 

Fig. 2. lES 1959+650 BeppoSAX spectrum of 28-29 
September 2001: LEGS + MEGS + PDS spectrum fit- 
ted with the curved model. The negative residuals left by 
simple and broken power-law models in PDS range disap- 
pear. 



by a constant model, with a maximum allowed variability 
lower than 20%. 

The second observation is longer (~ 5 x lO** sec) 
and can give us more information. During this run 
lES 1959+650 was still displaying a convex spectrum in 
the LEGS-MEGS range, with a break at - 2.5 keV and a 
further steepening in the PDS range. Therefore we again 
split our field of analysis extracting [0.15-2.5] keV LEGS 
and [2.5-10] keV MEGS fight curves, together with a PDS 
curve in the range [15-45] keV (Fig.OJ. For the LEGS and 
MEGS data we chose a binning time of 1000 seconds, while 
for the PDS we rebinned the data in 3600 second intervals, 
due to lower statistic. In each curve we accepted only bins 
with at least 20% of effective exposure time. Also dur- 
ing this observation BeppoSAX did not detect very large 
variability, with a possible indication of larger variability 
above the break. In fact, we get a high probability that 
the source is constant below the break for the LEGS data, 
while MEGS [2.5-10 keV] and PDS light curves are badly 
fitted by a constant model (constancy probability < 1%). 
This is not surprising since we are observing a steepening 
synchrotron spectra: a little hardening of the synchrotron 
component will cause small flux variations in the proxim- 
ity of the peak but larger flux variations as the spectrum 
steepens at higher energies. In any case these variations 
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Fig. 3. lES 1959+650 28-29 September 2001 fight curves. 
Top panel: [0.15-2.5] keV LEGS 1000 s binned light curve. 
Mid panel: [2.5-10] keV MEGS 1000 s fight curve. Bottom 
panel:[15-45] keV PDS 3600s binned light curve. 



are small and the main conclusion of this analysis is that 
the source did not clearly show brightness changes during 
the two BeppoSAX observations. Given that the source 
flux in the second observation increased by only ^^20% it 
is likely that it remained approximately stable at this high 
level at least from September 25 up to September 29. 



5. Optical observations 

We observed lES 1959+650 for about nine hours in the 
Cousins R band during the night of September 28-29 
2001 with the 50 cm f/4.5 Newtonian reflector of the 
Astronomical Station of Vallinfreda (Rome) (e.g. Maesano 
et al. 1997). The reference photometric sequence by 
Villata et al. (1998) was used, except for their star 5, 
which is variable as already pointed out by Nikolashvili, 
Kurtanidze and Richter (1999). The source remained con- 
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Fig. 4. The optical, Rc magnitude, light curve of lES 1959+650 since August 1995. Open triangles are data from the 
Perugia observations, filled circles are from Roma-Vallinfreda. The vertical line indicates the dates of the BeppoSAX 
pointings, activated by an X-Ray trigger. 



stant around a mean value of R=14.67 within our sensi- contribution of the underlying host galaxy. In the case of 
tivity of 0.02 mag (Maesano et al. 2002). 1ES1959+65 the host is detectable even with the short 

focal length of our 50 cm telescope. Heidt et al. (1999) 
To measure properly the optical SED of a blazar it is using the 2.5m NOT telescope at La Palma found that 
necessary to subtract from the observed flux values the 
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a De Vaucoulers profile (re = 9.5", Rtot=14.80) does not 
give an acceptable for the fit to the surface brightness 
of this galaxy, at variance with other host galaxies of their 
sample of BL Lacs. They found a better agreement using 
a different exponent (0.41 instead of 0.25), with an effec- 
tive radius re=10.3", a small additional disk component, 
with re=4.8" and a total galaxy R magnitude Riot=14.9. 
Furthermore they found a dust lane 1" north of the nu- 
cleus. Scarpa et al. (1999) using HST images confirmed 
the presence of the dust lane but found a satisfactory fit to 
the brightness profile with a De Vaucouleurs law (re=5.1", 
Rtot = 14.92) without the need for an additional disk. We 
remark however that the derived contribution of the host 
galaxy depends on the seeing conditions: a larger seeing 
spreads the galaxy image over a larger apparent radius, 
so that the galaxy appears fainter within a given radius. 
This effect is stronger if the intrinsic r,, is smaller than the 
seeing radius, which is not our case, however, whichever 
estimate of (4.8 or 9.5 arcsec) is assumed. 

For the purpose of this paper, i.e. to evaluate the 
contribution of the host galaxy within our photometric 
radius, we compute the R magnitude integrating a De 
Vaucouleurs profile for both estimates, finding R=15.66 
following Scarpa et al. (1999) and R=15.95 following Heidt 
et al. (1999) (see TableEJ. Both values are fainter than the 
minimum recorded luminosity of the source, as can be de- 
rived from the historical light curve shown in Fig. 0] built 
using the data obtained in the course of a monitoring pro- 
gram of a sample of bright blazars carried out since 1994 
by the Roma and Perugia groups. The telescopes used 
were the AIT (0.40 m) of the Perugia Observatory (Tosti 
et al. 1996) and the reflector telescope of Vallinfreda. 
The source shows a flux variability between R=14.4 and 
R=15.2, without evident periodicity. Some additional pho- 
tometric points have been published by Villata et al. 
(1998), also giving the source luminosity within the above 
quoted range. 

After subtraction of the galaxy contribution we com- 
puted the intrinsic flux of the source (in mJy) adopting 
the extinction curve by Ricke and Lebofsky (1985) and 
two possible values for the colour excess: E(B— V)=0.16 
(consistent with the Galactic value for the N^f, imply- 
ing A;j=0.43) and E(B-V)=0.26, consistent with the N// 
value derived from the X-ray data. The zero magnitude 
flux value (3.12 Jy at R=0.0) is taken from Elvis et al. 
(1994) (see TableE}. 

The V— R colour index of the source is remarkably 
constant in our database, 0.46, with a dispersion of 0.01. 
Thus, assuming an intrinsic value of B— V=0.61 for the 
host galaxy (observed B— V=0.77) following Fukugita et 
al. (1995), we can derive the V mag (see Table 13)) and a 
spectral energy slope value, in the Log(z/)-Log(Fiy) plane, 
of 0.49 and 0.83 for the two estimates of the host galaxy 
magnitude, respectively. 



6. The Spectral Energy Distribution 

We have constructed the SED of lES 1959+650 (shown in 
Fig. 5) using our iJeppoSAX and optical data (filled sym- 
bols) and data from the literature. We also include the pre- 
vious X-ray spectra measured by ROSAT and i?eppoSAX 
(from Beckmann et al 2002). The TeV spectrum measured 
during a large flare in May 2002 has been recently reported 
by Aharonian et al. (2003). The same paper gives the TeV 
flux observed for the period 2000-2001. In the SED we 
assumed that the TeV spectral shape remains almost un- 
changed with flux and simply rescaled the spectrum to the 
average flux observed during 2000-2001. 

The spectral data presented in the previous sec- 
tions show that the X-ray spectral distribution of 
lES 1959+650 is characterised by a well defined curva- 
ture. In the case of the parabolic model, a measure of this 
curvature is given by the parameter h in Eq.(2). It is inter- 
esting to compare this result with those of Mkn421, whose 
BeppoSAX. data have been analysed with the same spec- 
tral law (Massaro et al. 2003b). The typical b values for 
this source are higher than 0.3 and only during the bright 
phase of spring 2000 are close to 0.2, a value comparable 
to that measured in the September 28-29 observation of 
lES 1959+650 and equal to 0.21. If the spectral and lu- 
minosity trends of the sources are similar this may be an 
indication that in the previous observations, when it was 
much fainter (see Table 4), the spectral curvature could 
have been greater as also indicated by the steeper spectral 
slopes. Furthermore, with this model one can easily derive 
the energy of the maximum in the SED (see Massaro et 
al. 2003b for details) which is --0.1 and -0.7 keV in the 
two observations, respectively. Similar values are found in 
the intermediate luminosity states of Mkn421, while dur- 
ing flares the peak ranges between 2.5 and 5 keV. We 
can therefore conclude that the typical peak energies for 
lES 1959+650 are around ~0. 1-0.3 keV, approaching 1 
keV in the X-ray high states. 

We have used a homogeneous, one-zone synchrotron 
inverse Compton model to reproduce the SEDs of 
lES 1959+650. The model is very similar to the one de- 
scribed in detail in Spada et al. (2001), it is the "one-zone" 
version of it. The same model was applied to S5 0716+714 
and OQ530 by Tagliaferri et al. (2003); further details 
can be found in Ghisellini, Celotti & Costamante (2002), 
where the same model has been applied to low power BL 
Lacs. The main assumptions of the model are: 

— The source is cylindrical, of radius R and thickness 
AR' = R/T (in the comoving frame, where F is the 
bulk Lorentz factor). 

— The source is assumed to emit an intrinsic luminosity 
L' and to be observed at a viewing angle 6 with respect 
to the jet axis. 

— The particle distribution N{'^) is assumed to have the 
slope n [iV(7) oc 7~"] above the random Lorentz fac- 
tor 7c, for which the radiative (synchrotron and in- 
verse Compton) cooling time equals AR'/c. The mo- 
tivation behind this choice is the assumption that rel- 
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Table 5. Observed total (blazar+galaxy) R magnitude and derived intrinsic values for the blazar nucleus. The contri- 
bution of the galaxy is calculated assuming two different profiles (see text). The V magnitude has not been observed, 
but it is derived thanks to a very constant V-R value observed over the years for this source, see text. 



ativistic particles are injected into the emitting vol- 
ume for a finite time, which we take roughly equal 
to the light crossing time of the shell. This crossing 
time is roughly equal to the time needed for two shells 
to cross, if they have bulk Lorentz factor differing by 
a factor around two. The electron distribution is as- 
sumed to cut-off abruptly at 7max > 7c- We then as- 
sume that between some 7inin and 7c the particle dis- 
tribution A^(7) oc 7~("~i). This choice corresponds to 
the case in which the injected particle distribution is 
a power law (oc 7^("~i)) between ^rnin and 7niaxj with 
7min < 7c- Below 7min we assume ^^(7) oc 7"^. The 
value of 7max is not crucial, due to the fact that the 
electron distribution, for our sources, is steep. It has 
been chosen to be a factor ~ 100 larger than 7min- 

According to our X-ray data, the peak of the syn- 
chrotron spectrum occurs in the 0.1-0.7 keV band. This 
motivates the choice of the adopted f peak as listed in Table 
El This constrains the values of the input parameters of 
our model. Since we determine 7poak as the energy of 
those electrons that can cool in the injection time (i.e. 
7poak = 7c), there is a relation between z/peak and the 
cooling rate. If the latter is dominated by synchrotron 
emission, the location of t'peak constrains the value of the 
magnetic field. This in turn fixes the relative amount of 
inverse Compton radiation. 

The remaining degree of freedom for the choice of 
the input parameters is due to the value of the beam- 
ing factor (i.e. F and the viewing angle 9), and to the 
redshift. Another important input parameter is the size 
of the emitting region. Although the source did not show 
fast variability during our BeppoSAX observations, it is 
known to be variable, with a minimum observed vari- 
ability timescale of few hours (e.g. Holder et al. 2003a, b; 
Aharonian et al. 2003). The size of the emitting region is 
therefore constrained to be less than one light-day (i.e. 
R < ctvaid/i^ + z)). Although a one-zone homogeneous 
model is forced to use the minimum variability timescale 
observed at any band to constrain the size, it is also clear 
that this model is a simplification of a scenario which may 
be more complex. Our choice of i? = 9 x 10^^ cm and 
(5 '--^ 18 corresponds to a minimum variability timescale of 
tvar ~ 5 hours. 

The input parameters are listed in Table El and the 
model fits are shown in Fig.O The applied model is aimed 



at reproducing the spectrum originating in a limited part 
of the jet, thought to be responsible for most of the emis- 
sion. This region is necessarily compact, since it must ac- 
count for the fast variability shown by all blazars, espe- 
cially at high frequencies. The radio emission from this 
compact region is strongly self-absorbed, and thus the 
model cannot account for the observed radio flux. This 
explains why the radio data are above the model flts in 
the figure. 

7. Discussion and conclusions 

During both BeppoSAX. observations, triggered by an ac- 
tive X-ray status of the source, we detected X-ray spectra 
that steepen with increasing energies, up to 45 keV. The 
observed spectra are probably due to synchrotron emis- 
sion, with the synchrotron peak moving to higher energies 
with respect to previous observations. lES 1959-1-650 be- 
longs to the class of HBL, being characterized by a syn- 
chrotron peak in the soft X-ray range. This is confirmed 
by our BeppoSAX observations, which caught the source 
in a high X-ray state. The slope of the X-ray spectrum is 
harder than during previous X-ray observation (see Fig. 
0), suggesting a "harder when brighter" behavior, common 
to many other blazars, at least at frequencies above the 
synchrotron peak (indications exist that the same behav- 
ior is shown also at frequencies above the inverse Compton 
peak, see e.g. 3C279, Hallo et al. 2002; PKS 0528+134, 
GhiseUini et al. 1999). 

We do not know yet, however, if lES 1959+650, at 
still higher X-ray (2-10 keV) fluxes, is similar to the 
flaring state of Mkn501 (during summer 1997) with a 
flat (q!j; < 1) X-ray spectrum peaking at hundreds of 
keV, or resembles instead the X-ray spectrum of Mkn421 
which has been always observed with > 1. The phys- 
ical parameters derived by applying our one-zone SSC 
model are typical of all low power HHL objects, charac- 
terized by a relatively large beaming factor, low luminosity 
and absence of external seed photons testifled also by the 
very weak (if any) broad emission lines. One can compare 
the parameters derived in this paper with those derived 
by a sample of HL Lac objects in GhiseUini, Celotti & 
Costamante (2002) to see that all sources in this class are 
characterized by similar parameters. These are the sources 
which are among the best candidates to be strong TeV 
emitters (see Costamante & GhiseUini 2002), and in fact 
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Fig. 5. The SED of 1959+650. Non simultaneous data (star symbols) have been collected from the NASA Extragalagtic 
Database (NED). The "butterflies" in the X-ray region mark the ROSAT (lower) and the previous Beppo (upper) 
spectra (from Beckmann et al. 2000). The filled circles represent our optical and X-ray data. The lines refer to the 
SSC models for the two BeppoSAX observations. The TeV data (upper symbols) are from Aharonian et al. (2003) and 
refer to the flare in May 2002. For the plotted lower TeV data we assume no change in spectrum, and we have simply 
rescaled the May 2002 data by the amount corresponding to the average flux observed during the years 2000-2001, as 
reported by Aharonian et al. (2003). 



lES 1959+650 recently has been detected at TeV energies 
at levels above the Crab by the WHIPPLE and HEGRA 
Cherenkov telescopes (Holder et al. 2003a, Aharonian et 
al. 2003). At the time of the TeV flares the X-ray flux was 
larger than during our BeppoSAX. observations, in fact it 
was in the brightest X-ray state of the source in the last 
seven years, as monitored by ASM onboard RossiXTE 



(Holder et al. 2003b). This confirms the strong correla- 
tion between X-ray and TeV emission, in turn confirm- 
ing that the same population of electrons are emitting in 
both bands. However, one very rapid TeV flare (doubling 
timescale of 7 hours) was detected by WHIPPLE without 
a corresponding brightening of the X-ray flux (Holder et 
al. 2003b). This is clearly difficult to explain with one- 
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Table 6. Model Input Parameters. Column (1): name of the source; column (2): observation date; column (3): intrinsic 
luminosity i'; column (4): magnetic field B] column (5): size of the emitting region i?; column (6): bulk Lorcntz factor 
F; column (7): viewing angle 9 (in degrees); column (8): beaming factor d; column (9): slope of the particle distribution 
n; column (10): minimum Lorcntz factor of the injected electrons 7min; column (11): Lorentz factor of the electron 
emitting at the peaks, 7pcak; column (12): synchrotron peak frequency fpoak- Note that 7peak and :^pcak are derived 
quantities and not input parameters. 
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zone SSC model. If more examples of this intriguing flare 
episode will be detected in this and other sources, then 
some extra complexity in the model will be necessary. 

This is the second source in our BeppoSAX. ToO pro- 
gram that has been observed because it was in an active 
X-ray state, the other one being PKS2005-489 (Taghaferri 
et al. 2001). In both cases, only the synchrotron compo- 
nent was observed with the synchrotron emission peak- 
ing at 1-2 keV and no fast variability was detected. 
For six other sources the trigger came from the optical 
band. Thus, we have probably been biased toward sources 
that have an higher optical variability. These should be 
the blazars which have the synchrotron peak in the IR- 
optical band and these are of course essentially LBL or 
intermediate blazars. For all of them we detected either 
only the Compton component or both the synchrotron 
and Compton components (see Tagliaferri et al. 2002). 
Moreover in three cases we detected very fast variability 
in the X-ray band, but only for the synchrotron compo- 
nent (ON231: Taghaferri et al. 2000; BL Lac: Ravasio et al. 
2002; S5 0716-1-714: Taghaferri et al. 2003). This can be in- 
terpreted as the presence in the X-ray band of a Compton 
component (slowly variable on time scale of months) , and 
the tail of a synchrotron component with fast and the er- 
ratic variability. 

Thus, we can conclude that, although our sample of 
BeppoSAX ToO observations is limited, the behaviour of 
LBL-intermediate blazars in the X-ray band is different 
than that of HBL, when they are both in a high state of 
activity. 
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